Redistribution of particle and anti-particle entanglement in non-inertial frames 
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We analyse the entanglement tradeoff between particle and anti-particle modes of a Dirac field 
from the perspective of inertial and uniformly accelerated observers. Our results show that a redis- 
tribution of entanglement between particle and anti-particle modes plays a key role in the survival 
of fermionic field entanglement in the infinite acceleration limit. 
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I. INTRODUCTION 

Understanding entanglement in non-inertial frames has 
been central to the development of relativistic quantum 
information [l|-[l9| . The main aim of this field is to imple- 
ment quantum information tasks (such as quantum tele- 
portation) in relativistic settings. Quantum correlations 
are an important resource in most quantum information 
applications therefore, it has been of interest to under- 
stand how entanglement can be degraded [1, H, [l9l - l2ll | or 
created [22| - [25j by the presence of horizons or spacetime 
dynamics. 

It is well-known that the entanglement between modes 
of bosonic and fermionic fields is degraded from the per- 
spective of observers moving in uniform acceleration. In- 
terestingly, entanglement is completely degraded in the 
infinite acceleration limit in the bosonic case while for 
fermionic fields a finite amount of entanglement remains 
in the limit. However, the reasons for these differences 
were not completely clear. In this paper we show that a 
redistribution of entanglement between particle and anti- 
particle modes plays a key role for the preservation of 
fermionic field entanglement in the infinite acceleration 
limit. 

In our analysis we consider entangled states which in- 
volve particle and antiparticle field modes from the per- 
spective of inertial observers. Previous studies consid- 
ered entangled states involving exclusively particle modes 
from the inertial perspective. To study particle and an- 
tiparticle entangled states we develop a generalization of 
the formalism introduced in (26| which relates general 
Unruh and Rindler modes. This formalism refines the 
single-mode approximation [l|, [Tl| which has been ex- 
tensively used in the literature. In particular, we will 
consider in our analysis a fermionic maximally entangled 
state which has no neutral bosonic analog. This state 
which is entangled in the particle/ antiparticle degree of 
freedom can be produced, for example, in conjugated pair 
creation or in the production of cooper pairs. The analy- 
sis of such states is only possible under the mode transfor- 
mations we introduce here since the single approximation 
P, [ll| does not hold in this case. 

Considering a more general set of states from the in- 



ertial perspective allowed us to understand that in non- 
inertial frames entanglement redistributes between par- 
ticle and anti-particle modes. This is a somewhat sim- 
ilar effect to that observed in the inertial case: entan- 
glement redistributes between spin and position degrees 
of freedom from the perspective of different inertial ob- 
serves [27l . [28j . Interestingly, one can conclude that 
fermionic entanglement remains finite in the infinite ac- 
celeration limit due to this redistrubution of entangle- 
ment, which does not occur in the bosonic case. Our re- 
sults are in agreement with previous results which show 
that main differences in the behaivour of entanglement 
in the bosonic and fermionic case are due to Fermi-Dirac 
and Bose-Einstein statistics, contrary to the idea that he 
dimension of the Hilbert played an important role (29| . 

This paper is organized as follows: in section |ll] we in- 
troduce transformations between Minkowski, Unruh and 
Rindler modes for fermionic fields. This section extends 
results from[26| by including anti-particle modes. In sec- 
tion mil we analyse the entanglement transfer between 
the particle and antiparticle sectors in different kind of 
maximally entangled states when one of the observers is 
uniformly accelerated. Finally, conclusions and discus- 
sions are presented in section ITVl 



II. DIRAC FIELD STATES FORM THE 
PERSPECTIVE OF UNIFORMLY 
ACCELERATED OBSERVERS 

We consider a Dirac field in I + Idimensions. The 
field can be expressed from the perspective of inertial 
and uniformly accelerated observers. In this section we 
introduce the transformations which relate the mode op- 
erators and states from both perspectives. Such transfor- 
mations have been introduced in [2^ for particle states. 
Here we extend this results including transformations for 
anti-particle modes which will be needed in our analysis. 

Minkowski coordinates {t, x) are an appropriate choice 
of coordinates to express the field from the perspective for 
inertial observers. However, in the uniformly accelerated 
case Rindler coordinantes (77, x) rnust be employed. The 
coordinate transformation is given by 
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'Previously known as Fucntcs-Guridi and Fucntes-SchuUer. 



where < x < 00 and — 00 < 77 < 00. The transfor- 
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Figure 1: Rindler space-time diagram: lines of constant po- 
sition X = const, correspond to hyperbolae and all curves of 
constant 77 correspond straight lines that converge in the ori- 
gin. An uniformly accelerated observer Rob travels along a 
hyperbola constrained to either region I or region II. 



equation in Rindler coordinates) is given by 
i 

(4) 

where is again a normalisation constant. Cj^-^^dj^-^ 
with E = I, II represent Rindler particle/antiparticle op- 
erators. Note that operators in different regions E = I, II 
do not commute but anticommute. j = {Eq, s'} is again 
a multi-label including all the degrees of freedom. Here 
u^j and u^ii are the positive/negative frequency solu- 
tions of the Dirac equation in Rindler coordinates with 
respect to the Rindler timelike Killing vector field in re- 
gion I and II, respectively. The modes ("fc n) do not 
have support outside the right (left) Rindler wedge. The 
annihilation operators Cfe.M,c?fc.M define the Minkowski 
vacuum \0)-^ which must satisfy 

Ck,M \0)m = ^fe,M \0)m ^ 0' (5) 

In the same fashion Cj^^,dj,^, define the Rindler vacua 
in regions E = I, II 

c,,E |0)s = |0)s = 0, Vi, S = I, II. (6) 

The transformation between the Minkowski and Rindler 
modes is given by 



mation is defined in two spacetime regions \t\ < x and 
2: < — |i| called region I and II, respectively . The curve 
X = 1/a, where a is a positive constant of dimension in- 
verse length, is the world line of a uniformly-accelerated 
observer whose proper acceleration equals a. The proper 
time of this observer is given by rj/a in region I and by 
— r//a in region II. Note that djj is a timelike Killing vec- 
tor in both I and II, and it is future-pointing in I and 
past-pointing in II. 

The Dirac field satisfies the equation {i'y>^{d^ — F^) + 
m}0 = where 7'' are the Dirac-Pauli matrices and F^ 
are spinorial affine connections. The field expansion in 
terms of the Minkowski solutions of the Dirac equation 



(2) 



where Nm is a normalisation constant. The label ± 
denotes positive and negative energy solutions (par- 
ticles/antiparticles) with respect to the Minkowskian 
Killing vector field dt- The label k is a multilabel in- 
cluding energy and spin k ~ {E^^, s} where s is the com- 
ponent of the spin on the quantisation direction, and 
dk are the particle/antiparticle operators that satisfy the 
usual anticommutation rules 



{cfc,M,c|,' m} = {dk,Mi'^i'M} — ^kk', 



(3) 



and all other anticommutators vanishing. The Dirac field 
operator in terms of Rindler modes (solutions of the Dirac 



The coefficients which relate both set of modes are given 
by the inner product 



(ufe,u-,) = j d^xu\uj, (7) 
so that he Bogoliubov coefficients yield [1^ [13, HH , 



2V^ 



I* 



a" = («!)* /3" (/3I)* 7" = (7I)* 77" = {rfy 

(8) 

where tan ro = e , Eq is the energy of the Rindler 
mode fc, Eui is the energy of the Minkowski mode j and 9 
is a parameter defined such that it satisfies the condition 
En — m cosh 9 and jfeol — msinhfl (see [1^). Finally, 
taking into account that cj^m = {u'^M^'f') '^]m ~ 
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K,M''/')' find 



fe 



(9) 



We now consider the transformations between states in 
different basis. For this we define an arbitrary element 
of the Dirac field Fock basis for each mode as 



\F, 



where 



/II ' 



ip)r \i)u 



(10) 



(11) 



Here the ± signs denote particle/antiparticle. It is now 
convenient to introduce a new basis for inertial observers 
which corresponds to a superposition of Minkowski 
monocromatic modes. The reason for this is that the 
new modes, called Unruh modes [2^, and Rindlcr modes 
have a simple Bogoliubov transformation: each Unruh 
mode transforms to a single frequency Rindler mode. 
This transformation is given by 



Cfc3 = [cos Tk Cfc,i - sin „ 
Cfc,L = (cosrfc Cfcji - sinrfc d\j 
^k.K ^ (sinrfc Cfe.i + cosr^ 4ji 
^fe.L = ( si'^ '^k Ck.ii + cos Tk d\ -A , 



(12) 



were Ck^R,L and Dk^n^h are the Unruh mode operators. 

The corresponding transformation between Minkowski 
and Unruh modes are given by, 



k L 



where 



k L 



x(<r,®1l)+X*(1r®<l) 



X = il + i)e 



iOEn 



(13) 
(14) 

(15) 



Here wc have written the tensor product structure (|10p 
explicitly. 

For massless fields it can be shown [32| that the Unruh 
operators have the same form as Eq. (jl2p however in this 
case tanrfc = g-'^^/'^. 

In the massless case, to find the Minkowski vacuum in 
the Rindler basis we consider the following ansatz 



|o)m = (8) Io^^)m > 



(16) 



where |0o)m = |On)R ® |Oa)L. Wc find that 

|On)R, = E \nn,s)i \nn-s)u) 

n.s 



|0o)i 



E 

n.,s 



(Cn.n.s \nn,s)i |"-n,-s)n) , (17) 



where the label ± denotes particle/antiparticle modes 
and s labels the spin. The minus signs on the spin label in 
region II show explicitly that spin, as all the magnitudes 
which change under time reversal, is opposite in region I 
with respect to region II. 

We obtain the form of the coefficients Fn^n,s, Gn,n,s for 
the vacuum by imposing that the Minkowski vacuum is 
annihilated by the particle annihilator for all frequencies 
and values for the spin third component. 

Since the simplest case that preserves the fundamental 
Dirac characteristics corresponds to Grassnian scalars, 
we study them in what follows. In this case, the Pauli 
exclusion principle limits the sums (jl7p and only the two 
following terms contribute 

|0o)r, = ^i|0o)i+|0o)„+i^i|lo)i+|lo)„, 
|0o)l = Go|Oo)r|Oj2)A + Gi|lo)r|lo)iV (18) 

Due to the anticommutation relations we must introduce 
the following sign conventions 



|lf2)i^|lo)n 



4i4iil0o)i+|0o)„, 

\ln)i\ln)ti = dl,clji\On)^\On)ti, 

= -4xivdli\0n)i\0n)ij. (19) 

For the case of Grassman scalars [l^ , after imposing that 
Clu,m \On)M = lo'^ ^ omega we obtain Cq.r lOn)^ = 
Cn,L l^n)]^ = for all fl. The vacuum state then yields 



Oq) = [cosrn \On)i |Oo)ii + sinro \ln)i |1q)„ 



Using equation (|T4|) that this vacuum state also satisfied 
dui,M \On)M = which is equivalent to Dn.R |0o)r = 
Dn,i^ |On)L = Vfi. For convenience, wc introduce the 
following compact notation. 



\ijkl)n = \^n)t \jn)i |fco)^ |^a)„ 



(21) 



which is slightly different from the one employed in [261 . 
In this notation the vacuum state is written as, 

\Qn) = cos2ra|0000)j^-sinr,cosro|0110)j^ (22) 
+ sinrocosro |1001)jj - sin^ |llll)f2 . 

The Minkowskian one particle state is obtained by ap- 
plying the creation operator of particle or antiparticle to 
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the vacuum state = cj^ ^ |0)m, |1j7u 



|0), 



where the Unruh particle/antiparticle creator is a com- 
bination of the two Unruh operators 

4xv =PR.(^aR ® 11) +Pl(1r ® i??2,L)- (23) 
9Rj9l,Pr,Pl ai'c complex numbers satisfying |(7rP + 

kLp = l, |pRp + bLp = l. 

The parameters pr.l are not independent of (7r.l- We 
demand that the Unruh particle and antiparticle oper- 
ators are referred to particle and antiparticle modes in 
the same Rindler wedges. Therefore to be coherent with 
a particular election of and qh, we have to choose 
PL = qn and pn = qh, 



~-qR{Cl^<»lL)+qL{tR®Cl^), 



dlu =9L(^aR ® 1l) + gR(lR ® (24) 
The Unruh L and R field excitations are given by 



|lo)i |0n)„ 



\ln)t^Cl^\On)^^\On)i\ln)ti 
|1o)r =^sVr|0^2)r = |Oj^)+|lo)ri 

=<LlO^^)L = |lf2)r|Osi)+ (25) 

and therefore, 

=<ZR |lf2)R ® |0o)l + |0o)r ® , 

|lfe)u ='ZL |1o)r ® |0o)l + 'ZR |0f2)R ® . (26) 
In the short notation we have introduced the state reads, 

= gR[cosrfe|1000)j, -sinrn|1110)f,] 
+ qL [cosrn 10010) + sin ro |1011)f^] , 

|lfc)u = gL[cosrfc|0001)f,-sinro|0111)J 

+ gR[cosrn|0100)o+sinro|1101)jj]. (27) 



III. PARTICLE AND ANTI-PARTICLE 
ENTANGLEMENT IN NON-INERTIAL FRAMES 

Having the expressions for the vacuum and single par- 
ticle states in the Unruh and Rindler bases enables us to 
analyse the degradation of entanglement from the per- 
spective of observers in uniform acceleration. Let us con- 
sider the following maximally entangled states from the 
incrtial perspective 



* + ) =^(|0.)M|Oo)u + |l.)M|ln);), 



I*-) =^(|0..)M|Oo)u + |l-)M|li^>uj' 

|vl/i)=i=(|l^)+|lj,)^ + |l^)^|l,,) + 



(28) 
(29) 
(30) 



where the modes labeled with U are Grassman Unruh 
modes and the label cr = ± denotes particle and an- 
tiparticle modes. The first two states correspond to en- 
tangled states with particle and antiparticle Unruh ex- 
citations, respectively. These two states are analogous 
to the bosonic state --i=(|0)jyj^ \0)jj + |1)^^ \^)u) which is 
entangled in the occupation number degree of freedom. 
The third state has no analog in the neutral bosonic sce- 
nario since in this case the state is entangled in the parti- 
cle/antiparticle degree of freedom. In spite that fermionic 
entanglement in non-inertial frames has been extensively 
studied in the literature Q , states p9| and (pO)) have not 
been considered before. 

We consider Alice to be an inertial observer with a 
detector sensitive to w modes while her partner Rob 
who is in uniform acceleration carries a detector sen- 
sitive to n modes. To study the entanglement in the 
states from their perspective we must transform the 
modes to Rindler modes. Therefore, Unruh states must 
be transformed into the Rindler basis. The state in the 
Minkowski-Rinder basis becomes effectively a tri-partite 
system. As it is commonplace in the literature, we define 
the Alice- Rob bi-partition as the Minkowski and region I 
Rindler modes while the Alice- AntiRob bi-partitions cor- 
responds to Minkowski and region II Rindler modes. To 
study distillable entanglement we will employ the nega- 
tivity A/", defined as the sum of the negative eigenvalues 
of the partial transpose density matrix. Two cases of 
interest will be considered. In the first case we assume 
that Alice and Rob have detectors which do not distin- 
guish between particle and antiparticlcs. In this case, 
particles and antiparticlcs together are considered to be 
a subsystem. In the second case we consider that Rob 
and AntiRob have detectors which are only sensitive to 
particles (antiparticlcs) therefore, antiparticle (particle) 
states must be traced out. Our results will show that 
when Rob is accelerated, the entanglement redistributes 
between particles and antiparticlcs as a function of his 
acceleration. This effect is a unique feature of fermionic 
fields and plays an important role in the behavior of 
fermionic entanglement in the infinite acceleration limit. 



A. Entanglement in states ]'!'+) and 

To compute Alice-Rob partial density matrix in 
we trace over region II in |5'+)('I'+| and obtain. 



Par 



= ^[C^ |000)(000| + S'2C2(|010)(010| + |001)(001|) 

+ 5^ |011)(011| + |gRp(C2 |110)(110| + 52 |111)(111|) 

+ |gL|'(^'|110)(110| + C2|100)(100|) + g^(C3|000)(110| 
+S'2C|001)(111|) - g^(C2S'|001)(100| + 53|011)(110|) 
-gRg£5C| 111) (1001 ] +(H.c.)„^„_^.^g (31) 

were C — cos rn and S — sin . 
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The density matrix for the Ahce-AntiRob modes is ob- 
tained by tracing over region I, 



^4 



|000)(000| + S'2C2(|001)(001| + |010)(010|) 

|011)(011| + |gRp(C2 |100)(100| + S"^ |110)(110|) 

|(ZLp(52|lll)(lll|+C2|110)(110|)+<z^(C3|000)(110| 
S'2C|001)(111|) +qS,(C^'5'|001)(100| +53|011)(110|) 



(H.C.)non- 

diag. 



(32) 



To calculate the entanglement considering that Rob 
and AntiRob are able to detect both particles and an- 
tiparticles (the calculations follows from [1^), we first 
obtain the Alice-Rob and Alice- AntiRob partial trans- 
pose density matrices and their negative eigenvalues. The 
partial transpose matrix (pJifl)^^ is block diagonal and 
only the following two blocks contribute to negativity, 



Basis: {|100) , |010) , |111)} 

-qmlSC -qlS^ 
Basis: {|000) , |101) , |011)} 



(33) 



-qi^C^S 
-qlC^S q*^S^C 
q^S^C 



(34) 



Rob 



were the basis used is \iik) 



IM 



. Notice that 



although the system is bipartite, the dimension of the 
partial Hilbert space for Alice is smaller than the dimen- 
sion of the Hilbert space for Rob, which includes both 
particle and antiparticle modes. The eigenvalues only 
depend on IqrI and not on the relative phase between (/r 
and gL- 

We can carry out a similar calculation for the Alice- 
AntiRob subsystem. In this case we need to compute 
{p\fj)^'^ the blocks of the partial transpose density ma- 
trix which contribute to the negativity are, 

• Basis: {|111) , |010) , |100)} 



1 S\^. C^S^ q^C' |, 

2 1 



S'q*^ q*nqLSC\ 
qnqlSC qlC^ \qn\^C^ J 
Basis: {|011) , |101) , |000)} 



(35) 



54 qLS^C 
qlS^C q^^C^S 
qnC^S 



(36) 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 



Figure 2: Full Negativity for the Alice-Rob (Blue continuous) 
and Alice- AntiRob (red dashed) bipartitions as a function of 
rn = arctan e"'^^''" for various choices of |gR|. The blue con- 
tinuous (red dashed) curves from top to bottom (from bottom 
to top) correspond, to \q-R,\ — 1,0.9,0.8,0.7 respectively. 



Anti-Rob 

where we have considered the basis \ijk) = |fc)jj. 
Once more, the eigenvalues only depend on |gR| and not 
on the relative phase between (/r and q^- 

We see that, as discussed in [2^, when the accelera- 
tion increases the entanglement between Alice- AntiRob 
modes is created compensating the entanglement lost be- 
tween Alice-Rob. If < 1 the entanglement lost is 
not entirely compensated by the creation of entangle- 
ment between Alice- AntiRob resulting in a less entangled 
state in the infinite acceleration limit. An analysis of the 
quantum entanglement between Alice's modes and parti- 
clc/antiparticlc modes of Rob and AntiRob will be useful 
to disclose why correlations present this behaviour. 

We now analyse the entanglement when Rob's and 
AntiRob 's detectors are not able to detect antiparticlcs. 
In this case the entanglement is between their particle 
modes and Alice's subsystem. Since Rob cannot detect 
antiparticlcs we must trace over all antiparticle states 
and therefore, dSJ): = E„=o,i ("li" Par I">i ■ This 



yields 



Par+ 



1 



|00)(00| + |01)(01| +(7^^C|00)(11|) 



+ (kRl' + \<1l\'S') |11)(11| + |gL|'C^|10)(10|) 
+ (H.c.) .. (37) 

^ 'non-diag. ^ ' 

which is the partial state of Alice and the particles sector 
of Rob. 

The partial transpose (pJifl+)^'^ has only one block 
whose negative eigenvalue contributes to negativity 
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Basis: {|10) , |01)} 



1 f \qL\'C^ q*nC 

2 I QrC 



(38) 



The same procedure can be carried out for the system 
AR tracing over the antiparticle sector in ([5^ obtaining 



1 
2 



|00)(00| + 5'^ |01)(01| + g^C|00)(ll| 

qn\'s')\n){n\ + \qn\'c'\w){io\ 

(39) 



(H.C.)non- 

diag 



The partial transpose (p^;j+ Y'^ has only one block whose 
negative eigenvalue contributes to negativity 



Basis: {|10) , |01)} 



1 { \qn\^C^ qlC 

2 I qLC 52 



(40) 




therefore, 

+ (kLp + |<7RpC2) |10)(10| + \qn\'S' |11)(11|) 



|00)(00| + |01)(01| - qlS |01)(10| 



(H.c. 



non-diag. 



(41) 



is the partial state of Alice and the particles sector of 
Rob. 

The only block giving negative eigenvalues is 



Basis: {|11) , |00)} 



1 f \qn\^S^ ~qlS 

2 I -qLS 



(42) 



The density matrix for the Alice- AntiRob antiparticle 
modes is obtained by tracing over the particle sector in 



Par- 



1 r 

2 



C2|00)(00| + S'2|01)(01| 

(|gRp + kpc2)|io)(io| + |gLp^'|ii)(ii| 

(43) 



+ g^5|01)(10|) +(H.c.)„o„- 

J diag. 

Again only one block of the density matrix contributes 
to negativity 

• Basis: {|11) , |00)} 



q^S 



(44) 



The analysis of entanglement in the state (|29p is done 
in a completely analogous way. We find that the entan- 
glement behaves exactly the same way as in state ((28)) 
only that the role of particles is replaced by anti-particles. 
Therefore negativities are related in the following way 



^Ir^^^ar- ^1r. 



AR- 



AR+ 



Kr- 



■^1r- 



■^AR^ 



Figure 3: Negativity in the particle sector (Rob and An- 
tiRob can only detect particle modes) for the bipartition 
Alice-Rob (Blue continuous) and Alice- AntiRob (red dashed) 
for the state (|28|) as a function of rn = arctan e"'^"'^" for 
various choices of \qR\. The blue continuous (red dashed) 
curves from top to bottom (from bottom to top) correspond, 
to I^rI = 1,0.9,0.8,0.71 respectively. For \qR\ = 1 Alice- 
AntiRob curve is zero Va 

A similar calculation can be carried out considering 
that Rob and AntiRob detectors are only sensitive to 
antiparticles, i.e. tracing over particle states. In this 



case we obtain, ([HH): pj^^- = J2n=< 



1+ 



0,1 



Par I 



and 



^1r 



AR 



Kr = ^'ar 



We sec in figure [5] that the total entanglement between 
Alice and AntiRob starts decreasing and presents a min- 
imum before starting to grow again for higher accelera- 
tions. If \qji\ < 1 the entanglement in the limit a — ^ 
is distributed between the bipartitions AR and AR . 
The entanglement lost in the bipartition AR is not en- 
tirely compensated by the creation of entanglement in 
AR and therefore, this results in a state containing less 
entanglement in the infinite acceleration limit. 

Interestingly, the correlations between Alice and the 
particle sector of Rob and AntiRob always decrease (Fig. 
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0.25 r 




Figure 4: Negativity in the antiparticle sector (Rob and An- 
tiRob can only detect antiparticle modes) for the bipartition 
Alice-Rob (Blue continuous) and Alice- AntiRob (red dashed) 
for the state (|28|) as a function of rn = arctan e^'^^''" for 
various choices of |gR|. The red dashed (blue continuous) 
curves from top to bottom (from bottom to top) correspond 
to \q-R,\ = 1,0.9,0.8,0.71 respectively. 



[3|) while the correlations between Alice and the antipar- 
ticle sector of Rob and AntiRob always grow (Fig. H]). 
This behaviour explains why entanglement always sur- 
vive the infinite acceleration limit for any election of qi^ 
and q^. As Rob accelerates there is a process of entangle- 
ment transfer between the particle and antiparticle sector 
of his Hilbert space. The same happens with AntiRob, 
such that neither for AR nor AR the entanglement van- 
ishes for any value of the acceleration. 

For the simplest case \qj^\ = 1 we see that all the en- 
tanglement is initially (a — > 0) in the particle sector of 
the bipartition AR. As the acceleration increases the en- 
tanglement is transferred to the antiparticle sector of the 
bipartition AR such that, in the limit of infinite acceler- 
ation entanglement has been equally distributed between 
these two bipartitions. 

The tensor product structure of the particle and an- 
tiparticle sectors (|20p plays an important role in the be- 
havior of entanglement in the infinite acceleration limit. 
In the case of neutral scalar fields there are no antipar- 
ticles and entanglement is completely degraded. Note 
that in the case of charged bosonic fields there arc in- 
deed charged conjugate antiparticles. However, in this 
case the Hilbert space has a similar structure to the un- 
charged field (3]. The existence of bosonic antiparticles 
simply adds another copy of the same Hilbert space and 
no entanglement transfer is possible between particle and 
antiparticles. Clearly, the Hilbert space structure in the 
fcrmionic case (|20p is different and thus, we observe the 
differences in the entanglement behavior. 



Now, if we move to less trivial cases where Qr 1 the 
situation gets more complicated. In these scenarios we 
initially start with some entanglement in the particle sec- 
tor of AR, and there can also be an entanglement transfer 
between this sector and AR antiparticle sector. However 
no entanglement at all is transferred to the antiparti- 
cle sector of the subsystem AR unless the acceleration 
reaches a threshold given by 



2 m 

cos r = - — - 



1. (45) 



The maximum value of cos [r{a) 

|2 



is cos [r(a — > oo)] 



1/a/2 therefore for \qR\'^ > 2/3 entanglement is not trans- 
ferred to the antiparticle sector of Alice- Rob for any value 
of the acceleration. This explains why when q^ ^ 1 
the entanglement loss in the bipartition AR in the limit 
a — oo is smaller than in the extreme case qn = 1. 

It is evident that the choice of Unruh modes influence 
the transfer of entanglement between particle and an- 
tiparticle sectors. When the acceleration is larger than 
(PSj). when Igpil grows closer to 1/V2 more entanglement 
is transferred from the particle sector of AR to the an- 
tiparticles of AR. In the limit = 1/a/2 the same 
amount of entanglement is transferred to the antiparticle 
sector of both AR and AR. 

The particle entanglement (always monotonically de- 
creasing) resembles the behaviour of bosonic entangle- 
ment studied in (26| . Bosonic entanglement is monotoni- 
cally decreasing for both AR and AR subsystems. In the 
bosonic case there are not antiparticles and, hence, there 
is no possibility of entanglement transfer to antiparticle 
sectors of AR. This is the origin of the differences in en- 
tanglement behaviour for neutral bosons and fermions. 



B. Entanglement in state |3'i) 

This state has no neutral bosonic analog since it is 
entangled in the particle/antiparticle degree of freedom. 
Therefore, the analysis of entanglement in this state rev- 
els interesting features which are of genuinely fermionic 
nature. 

To study this type of state we em ploy our general- 
ization of the formalism developed in [2^] which relates 
general Unruh modes with Rindler modes. This formal- 
ism refines the single-mode approximation introduced in 
[J [ll| which has been extensively used in the literature. 
For this type of state the single approximation used in Q 
does not hold and attempting to use it leads to misleading 
results: one finds that maximally entangled states from 
the incrtial perspective appear disentangled from the ac- 
celerated perspective, irrespectively of the value of ac- 
celeration. Using the mode transformation introduced in 
our first section leads to sensible results: acceleration be- 
haves regularly for accelerated observers and approaches 
a maximally entangled state in the inertial limit. 

For convenience we will introduce a new notation for 
this case. For Alice, we denote the states by |-|-) if they 
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correspond to particles and |— ) for antiparticles. There- 
fore the state is written as 

I*i> = ^(|+>m|1o)u + |->m|1^^)J) (46) 

The density matrix for the subsystem Ahce-Rob is ob- 
tained from tracing over region II 

p\r -kRl'C^ l + 10>(+10| + \qK?S^ | + 11>(+11| 

+ |gLpC2 l+OOX+OOl + I^Lp^' l+10)(+10| 
-SCqnql 1 + 11) (+00| + |<7LpC2 |_00)(-00| 
+ |gLp5' h01)(-01| + IgRpC^ h01)(-01| 
+ \qK?S^ |-11)(-11| + SCq^qi^ |-00)(-ll| 
+ (|gRpC2 - I^Lp^") |+10)(-01| + (H.c. 



non- 
diag. 




% 0.1 0.2 0.3 0.4 0.5 0.6 0.7 



(47) 



and for the AHce-AntiRob partition we obtain 

p\r ^kRp^^' 1+00) (+00| + |gRp^2 1+10) (+10| 
+ |gLpC2 |+10)(+10| + |gLp52 |+11)(+11| 
+SCqKql |+00)(+ll| + |<7LpC2 |_01)(-01| 
+ |<?l|'52 |-11)(-11| + \q^\^C^ |-00)(-00| 
+ |gRp52 |-01)(-01| - SCqlqi^ hll)(-00| 
+(|gLrc2 - IgRp^^) |+10)(-01| + (H.c. 



non- 
diag. 

(48) 

Assuming that the observers cannot distinguish be- 
tween particle and antiparticles yield matrices where only 
one block of the partial transpose density matrix gives 
negative eigenvalues 

• Basis; {|-10) , |+01)} 



1 f (|(7RpC2 - \q^\^S^) 

(l9RpC2-|gL|252) 

In this case the negativity is given by 



(49) 



A similar result is obtained for the system Alice- 
AntiRob p\pc In this case the only block of the partial 
transpose that contributes to negativity is 

• Basis: {|-10) , |+01)} 



Figure 5: Particle-antiparticle maximally entangled state 
(|30p : Negativity for the bipartition Alice-Rob (Blue con- 
tinuous) and Alice- AntiRob (red dashed) as a function of 
rn ~ arctan e^'^^''" for various choices of |gR|. The blue con- 
tinuous (red dashed) curves from top to bottom (from bot- 
tom to top) correspond, to = 1,0.9,0.8,0.7 respectively. 
The curves for Alice- AntiRob entanglement have a minimum 
where the negativity vanishes. 



Interestingly, when Rob and AntiRob arc not able to 
detect either particle or antiparticle modes the entangle- 
ment in the state vanishes. The partial density matrices 
for AR and AR in this case yield 



Par+ 



=(kRp 

-(H.c.) 



|gLp52)|+l)(+l|- 
|gRpC2)|-0)(-0| 



non- 
diae. 



kLpC2|+0)(+0| 
|gRp^2|-l)(-l| 
(51) 



p\r- =(I9LP + IteP^') | + 1)(+1| + \qK?C^ l+0)(+0| 

+(|gRp + kLpc^) |-o)(-o| + kLp^' |-i)(-i| 

+ (H.c.)„„„. (52) 



p\r+ 



-(H.c.) 



kRpC2)|+0)(+0| 
\q^\^S')\-l){-l\- 



non- 
diae;. 



kpc2|-o)(-o| 

|gRp52| + l)(+l| 
(53) 



1 

2 UI9LPC2-|9R|252) 

resulting in 



qi.?C^-\qK?S^) , 



p\r+ 



-o)(+o| + kRrc^|-o)(-o| 



= (|gR|' + \qi.?C^ 

+ (I'ZLP + |<ZRp5') |-1)(-1| + \qi.?S^ | + 1)(+ 



|9L 

-(H.c.) 



non- 
diag. 



i)(+i| 

(54) 



for which negativity is strictly zero. The entanglement 
in this state is of a different nature as the entanglement 
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in states |^'+) and l^*-) therefore, a direct comparison of 
the behaviour of entanglement cannot be done. The total 
entanglement here is associated to correlations between 
particles and antiparticles and, therefore, if we trace out 
either the particles or antiparticles we effectively remove 
all the correlations codified in this degree of freedom. 

In the case when the detectors do not distinguish be- 
tween particles and antiparticles (see Fig. [5]) we found 
that the entanglement in the Alice- AntiRob bipartition is 
degraded with acceleration vanishing at a critical point. 
For higher accelerations entanglement then begins to 
grow again. Namely, the entanglement on the biparti- 
tion Alice-AntiRob, vanishes for a specific value of the 
acceleration if | | < 1 . This value of the acceleration is 
given by 



What is more, the surviving entanglement in the limit 
a — >■ oo is 

Therefore, when Igpl = I^lI = l/\/2 no entanglement 
survives in the limit of infinite acceleration. 

This shows that entanglement has a non-vanishing 
minimum value in the infinite acceleration limit (regard- 
less the election of Unruh modes) only when there is 
transfer of entanglement between particles and antipar- 
ticles. Otherwise, it is possible to find an Unruh mode 
whose entanglement vanishes in the infinite acceleration 
limit as in the bosonic case. We therefore conclude that 
the entanglement transfer between particle and antipar- 
ticle sectors plays a key role in explaining the behaviour 
of entanglement in the infinite acceleration limit. 

IV. CONCLUSIONS 

Including antiparticles in the study of fermionic en- 
tanglement allowed us to understand key features which 



explain the difference in behaviour of entanglement in 
the fermionic and bosonic case. Namely, we have shown 
that there is an entanglement redistribution between the 
particle and antiparticle sectors when Rob is in uniform 
acceleration. This entanglement transfer is not possible 
in the bosonic case and, therefore the differences in the 
behaviour of entanglement in the bosonic and fermionic 
case arise. In particular, we have shown that this entan- 
glement tradeoff gives rise to a non-vanishing minimum 
value of fermionic entanglement in the infinite accelera- 
tion limit for any choice of Unruh modes. 

We also exhibit a special fermionic state for which 
entanglement transfer between particle and antiparticle 
states is not possible. Interestingly, in this case we can 
find a specific choice of Unruh modes such that entan- 
glement vanishes in the infinite acceleration limit. In- 
cidentally, this choice (|gFi| = I^lI = l/"v/2) minimises 
the surviving entanglement of states and (P^]) . We 
showed that it is the tradeoff between the particles and 
antiparticles sector what protected them from a complete 
entanglement loss. 

Our analysis is based on an extension to antiparticles 
of the formalism introduced in [26j which relates Unruh 
and Rindler modes. This allowed us to analyse a more 
general family of fermionic maximally entangled states 
for which the single-mode approximation does not hold. 

This study sheds light in the understanding of rel- 
ativistic entanglement: the differences in bosonic and 
fermionic statistics give rise to differences in entangle- 
ment behaviour. This provides a deep insight on the 
mechanisms which makes fermionic entanglement more 
resilient to Unruh-Hawking radiation. 
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